    %

clear

diary('./table_5.txt');

for tau=25:25:75,

    filename='./data_results_tau';
    load([filename num2str(tau) '.mat'],'DWL*','x2_interest_univ');
    
        %%%%%%%%%%%%%%%%
        % DWL ESTIMATES
        %%%%%%%%%%%%%%%%

    fprintf('\n');
    fprintf('-------------------------------------------------------------------------------------------\n');
    fprintf('tau=%4.2f        unconstrained              constrained                linear quantile\n',tau/100);
    fprintf('-------------------------------------------------------------------------------------------\n');
    fprintf('  income        DWL  DWL/tax  DWL/inc      DWL  DWL/tax  DWL/inc      DWL  DWL/tax  DWL/inc\n');
    fprintf('-------------------------------------------------------------------------------------------\n');

    for yy=3:-1:1,

        fprintf('%8.0f   %8.2f %8.2f %8.2f %8.2f %8.2f %8.2f %8.2f %8.2f %8.2f \n', ...
		    round(exp(x2_interest_univ(yy,1))),DWL_unconstr(yy,1),10^2*DWL_pertax_unconstr(yy,1),10^4*DWL_perinc_unconstr(yy,1), ...
            DWL_constr(yy,1),10^2*DWL_pertax_constr(yy,1),10^4*DWL_perinc_constr(yy,1), ...
            DWL_qlin(yy,1),10^2*DWL_pertax_qlin(yy,1),10^4*DWL_perinc_qlin(yy,1) ...
			);

    end;

    fprintf('-------------------------------------------------------------------------------------------\n');

    clear DWL* x2_interest_univ

end;

diary off 

    %
